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Abstract; [Aim] Eysarcoris aeneus, a main pest of many economic crops, is widely distributed in Palaearctic 
region with humeri variable. Difference in the humerus length of “rounded humerus form" and “needle-like 
humerus form" suggests that intraspecific variation may exist in E. aeneus. [Methods] In this study, 142 
specimens (98 for E. aeneus and 44 for the outgroup E. ventralis ) from 19 regions of China were examined to 
ascertain the morphological variation of four characters (forewing, hindwing, head and scutellum) between the 
two forms, "rounded humerus form" and " needle-like humerus form" , respectively, of E. aeneus with the use 
of geometric morphometrics approaches, including multivariate regression, principal component analysis 
(PCA), canonical variate analysis (CVA) and discriminant function analysis (DFA). [Results] Allometry 
was not found in all of the four characters (forewing, hindwing, head and scutellum) examined. For the four 
characters, though the principal component analysis showed overlap between the two forms of E. aeneus, the 
canonical variate analysis showed significant differences between them ( P-values for Procrustes distance and 
Mahalanobis distance were all less than 0. 01). The correctly assigned percentages of “rounded humerus form” 
and “needle-like humerus form” ranged from 67% to 89% specimens based on the four characters using the 
discriminant function analysis. [Conclusion] This study proved that forewing, hindwing, head and scutellum 
can serve efficient evaluation criteria in describing intraspecific morphological difference of E. aeneus, and 
scutellum is the most informative character. The geometric morphometric descriptors showed significant shape 
differences and overlapped centroid size between the two forms of E. aeneus, indicating that the analysis of 
shape is more sensitive for detecting differences among lineages than that of size. Our study shows that 
geometric morphometrics is efficient in describing intraspecific morphological difference of E. aeneus, and lays 
a foundation for intraspecific variation studies of stink bugs in future. 

Key words: Eysarcoris aeneus; geometric morphometrics; principal component analysis (PCA) ; canonical 
variate analysis (CVA) ; discriminant function analysis (DFA) 
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variable. The “ needle-like humerus form’ 


INTRODUCTION 


doi; 10. 16380/j. kexb. 2019. 09. 009 


Eysarcoris aeneus, belongs to Eysarcoris 
( Hemiptera: Pentatomidae: Pentatominae ), is a 
pest of many economic crops such as wheat, 
soybean, rice, etc. (Ito, 1978; Learmonth, 1980; 
1985; Nasiruddin and Roy, 2012 ). 


According to some reports, E. aeneus, E. ventralis 


Zhang, 


and E. trimaculatus are main pests of rice, and the 
rate of pecky rice caused by E. aeneus at milk ripe 
stage was 10.2% (Lee et al., 2009; Nasiruddin and 
Roy, 2012). E. 


distributed in the Palaearctic region, with humeri 


aeneus is a species widely 
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described as a distinct species ( E. parvus) by Uhler 
(1896) , and then it was treated as a junior synonym 
of E. aeneus by Josifov and Kerzhner ( 1978). 
Difference in the humerus length of “ needle-like 


, 


humerus form” and ^" rounded humerus form’ 
suggests that intraspecific variation may exist in f. 
aeneus. 

Geometric morphometrics is a widely used 
statistical approach in morphology based on Cartesian 
landmark coordinates (Gilchrist et al., 2000; Debat 
et al., 2003; Mitteroecker and Gunz, 2009; Bai, 


2017). With the use of powerful and comprehensive 
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statistical procedures, geometric morphometrics can 
detect tiny shape differences of a morphological 
feature ( Rohlf and Marcus, 1993; Marcus et al., 
2000; Adam et al., 2004). Geometric morphometric 
analysis has proved 
discriminating sibling species, genus delimitation 


already been useful in 
and even in identification at the individual level 
among insects ( Matias et al., 2001; Baylac et al., 
2003; Dujardin et al., 2003; Villemant et al., 
2007 ; Sadeghi et al., 2009; Bai et al., 2014; Su et 
al., 2015; Giordani et al., 2017). 

Eysarcoris can be easily distinguished from 
other genera because of its basal angle feature in the 
scutellum with small, yellow, or pale smooth spots. 
Nevertheless, Eysarcoris species show complex 
variances ( Yang, 1962). In our former studies, 
geometric morphometric approach was employed to 
distinguish sibling species and the whole genus of 


Eysarcoris in China (Li et al., 2017a, 2017b). The 


aim of our present study is to quantify the conspecific 
morphological features between the two forms of E. 
aeneus based on the four characters ( forewing, 
hindwing, head and scutellum ) , and to test the 
possible use of these characters for intraspecific 
studies of Eysarcoris. 


2 MATERIALS AND METHODS 


2.1 Insects 
A total of 142 specimens (41 for " needle-like 
and 57 for 


form" ) of E. aeneus and E. ventralis were examined 


humerus form ” " rounded humerus 
in this study. The number and localities of the 
specimens were provided in Table 1. E. ventralis 
was chosen as the outgroup (Li et al., 2017b). 
Samples were collected mainly by manually pick-up 
or by using sweep nets from fields in China. All 
deposited 
University and Nankai University in China. 


specimens were in Taiyuan Normal 


Table 1 Collecting localities and body parts of Eysarcoris specimens used in geometric morphometric analysis 











Species Form Collecting locality Forewing Hindwing Head Scutellum 
Needle-like humerus form — Linyi, Shandong * 12913 8 14 913 ô 209218 20220 ô 

Pangquangou Nature Reserve, Shanxi " 8958 1897 ô 5958 4926 
Luya Mountain Nature Reserve, Shanxi * 16 8958 3928 2928 

"2 Qingfeng Mountain, Shaanxi * 19 29 

Rounded humerus form Taibai National Forest Reserve, Shaanxi 39 39 39 39 

Xunyang, Shaanxi 1918 19168 1918 1916 
Dongyang, Zhejiang 3918 39 
Huanan, Heilongjiang 4948 4948 4948 4948 

E. ventralis Beijing * 19 19 19 
Xintai, Shandong 39 29 39 39 
Linyi, Shandong * 10916 79 9928 99168 
Luo Mountain, Zhaoyuan, Shandong * 18 18 18 
Tai'an, Shandong * 16 18 18 
Xianren Island, Shandong * 39138 2938 39138 49138 
Huhuan Forest Park, Changsha, Hunan 19 19 19 
Cangyan Mountain, Hebei * 1926 1926 1928 
Yixian, Hebei * 18 16 1 ô 
Hongtong, Shanxi 16 16 16 16 
Baxian Mountain, Tianjin * 19 19 19 
Dongyang, Zhejiang 38 16 16 38 
Houhe Nature Reserve, Hubei * 19 19 19 


* The data of specimens of these localities are from Li et al. (2017a, 2017b). 


2.2 Geometric morphometric approach 

For all of the four characters ( forewing, 
hindwing, head and scutellum) , the same landmark 
selection method was used as the former works ( Li et 
al., 2017a, 2017b) ( Fig. 1). Landmarks and 
curves were based on homologous or corresponding 
criteria. Tps-DIG Ver. 2. 05 (Rohlf, 2006) was 


used to draw curves and landmarks and all semi- 


landmarks were converted to landmarks. All pictures 
and digital procedures were operated by the same 
person. 

Correlations between the Procrustes and tangent 
shape distances were calculated using tps-SMALL 
Ver. 1.20 (Rohlf, 2003). Centroid size (CS) was 
compared to evaluate the overall size variation of the 
four characters among groups. It was the square root 
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Fig. 1 Landmarks and curves selected 


> 


of the summed squared distance of each landmark 
from the center of the form ( Zelditch et al., 2004). 
A one-way ANOVA and HSD pairwise comparisons 
were conducted using paleontological statistics 
software Ver. 3. 06 ( PAST) ( Hammer et al., 
2001). To explore how shape varies with size, a 
multivariate regression of Procrustes coordinates 
against centroid size was conducted using MorphoJ 
1. 06d ( Klingenberg, 2011). The significance of 
the allometry was calculated with a permutation test 
with 10 000 iterations. 

To examine the shape variation of the four 
characters, the raw landmark data were firstly 
superimposed using the Procrustes methods in 
MorphoJ 1. 06d ( Bookstein, 1991; Klingenberge, 
2011). 


principal component analysis. The relative similarity 


Shape variations were assessed using a 


of the groups was analyzed using canonical variate 
( CVA ). 


Mahalanobis distance were calculated during the 


analysis Procrustes distance and 
CVA analysis. Discriminant function analysis with 
cross-validation was performed on the Procrustes 


superimposed coordinates. 


3 RESULTS 


Variation of the samples in shape space was 
strongly correlated with tangent space for all the four 
anatomical characters. For forewing, hindwing, 
head and scutellum, the correlation between the 


tangent space regressed onto Procrustes distance 


were 0. 999265, 0. 999623, 0. 999528 and 


: Forewing; B: Hindwing; C: Head; D; Scutellum. 


0. 999667 , respectively. 

For forewing and hindwing, the “ rounded 
humerus form” had a bigger centroid size than the 
“needle-like humerus form”, while for head and 
scutellum the centroid size of “ needle-like humerus 
form” was bigger than that of the “rounded humerus 
form” (Fig. 2). For all of the four characters, the 
difference in centroid size was nonsignificant between 
the two forms of E. aeneus ( Table 2). Allometry 
was not found in all the four characters of the groups 
since each multivariate regression was statistically 
nonsignificant ( P >0.05). 
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Fig. 2 Comparisons of centroid size of “rounded humerus 
form” and “needle-like humerus form” of 
Eysarcoris aeneus and E. ventralis 
Data in the figure are mean + SD. 
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Table 2 Tukey HSD [P-values (above) and studentized range statistic q (below) ] for 
the centroid size between the “rounded humerus form” and “needle-like humerus form” 
of Eysarcoris aeneus and E. ventralis 
Rounded humerus form Needle-like humerus form E. ventralis 
Forewing 
Rounded humerus form 0.5828 0.7537 
Needle-like humerus form 1.4050 0.9591 
E. ventralis 1.0150 0.3899 
Hindwing 
Rounded humerus form 0.2676 0.1493 
Needle-like humerus form 2.2070 0.9448 
E. ventralis 2.6610 0.4544 
Head 
Rounded humerus form 0.4827 0.0168 
Needle-like humerus form 1.6280 0.0003 
E. ventralis 3.8790 5.5070 
Scutellum 
Rounded humerus form 0.1253 0.0411 
Needle-like humerus form 2.1560 0. 8848 
E. ventralis 3.4220 0.6664 
The shape changes of the different groups were forewing, hindwing, head and scutellum, 
implied by the PCA scatter plots. The first two PCs respectively. The principal component analysis 


were plotted to indicate variation along the two axes 
and accounted for 73. 386% , 51. 301% , 82. 333% 
and 75. 017% of the variation among groups for 
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Principal component 1 


reveals overlap between the “ rounded humerus 
form” and “needle-like humerus form” specimens in 


all the characters (Fig. 3). 
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Fig. 3 Principal component analysis (PCA) of “rounded humerus form” and “needle-like humerus form" 
of Eysarcoris aeneus and E. ventralis 
A; Forewing; B: Hindwing; C; Head; D; Scutellum. 
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The CVA scatter plots of shape differences all groups were significantly different ( P < 0. 0001), 


showed that each group occupied different area and Procrustes distances were similar ( P «0.01) in 
(Fig. 4). Mahalanobis distances among the three all pairwise comparisons ( Table 3). 
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Fig. 4 Canonical variates analysis (CVA) of “rounded humerus form" and “needle-like humerus form" 
of Eysarcoris aeneus and E. ventralis 
A; Forewing; B; Hindwing; C; Head; D; Scutellum. 


Table3 Difference in Mahalanobis distance and Procrustes distance [ P-values ( above) , and distance ( below) ] between 
the “rounded humerus form” and “needle-like humerus form” of Eysarcoris aeneus and E. ventralis 
in the four characters (forewing, hindwing, head and scutellum) 

















Mahalanobis distance Procrustes distance 
Rounded humerus — Needle-like humerus : Rounded humerus — Needle-like humerus ; 
form form E. ventralis form form E. ventralis 

Forewing 
Rounded humerus form <0. 0001 <0. 0001 0.0016 <0. 0001 
Needle-like humerus form 18. 1086 <0. 0001 0.0182 <0. 0001 
E. ventralis 44. 6643 48. 4641 0.0574 0.0448 
Hindwing 
Rounded humerus form <0. 0001 <0. 0001 0.0010 <0. 0001 
Needle-like humerus form 16.9728 <0. 0001 0.0135 <0. 0001 
E. ventralis 39.2818 34.5077 0.0341 0.0312 
Head 
Rounded humerus form <0. 0001 <0. 0001 <0. 0001 <0. 0001 
Needle-like humerus form 5.8227 <0. 0001 0.0343 «0.0001 
E. ventralis 13. 7213 14. 3363 0.0627 0.0612 
Scutellum 
Rounded humerus form «0. 0001 «0.0001 «0.0001 «0.0001 
Needle-like humerus form 8.9854 «0.0001 0.0359 «0.0001 


E. ventralis 16.2703 18.1803 0.0233 0.0532 
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Almost all 
“ needle-like 


correctly assigned to their own groups, 


" rounded humerus form” and 
humerus form” specimens were 
when 


compared to E. ventralis. However, the percentages 


correctly assigned of “rounded humerus form” and 
“needle-like humerus form” specimens from cross- 
validation were mostly lower than 80% (Table 4). 


Table 4 Mahalanobis distances and correct grouping percentage in discriminant function analysis (DFA) 
( permutation test =1 000) of “rounded humerus form" (ER) and “needle-like humerus form” (EL) 
of Eysarcoris aeneus and E. ventralis (EV) 


Mahalanobis distance 


Percentage correctly assigned 











Compared to EV ER and EL Compared to EV 
ER and EL 
ER EL ER EL ER EL 
Forewing 7.1493 18.2040 15. 3005 6896 6796 94% 96% 
Hindwing 8.2789 14. 9685 10.4564 7696 7096 95% 96% 
Head 6.5737 22.4490 20. 3487 7196 7396 97% 100% 
Scutellum 10. 9763 12.1191 13.7435 89% 80% 89% 88% 


4 DISCUSSION 


Geometric morphometrics is sensitive to small 
variations and can detect subtle difference within 


CVA. 


discriminating populations, 


It was proved to be very powerful in 
groups and even 
subpopulation within the same location ( Mendes et 
al., 2007; Gidaszewski et al., 2009; Oleksa and 
Tofiski, 2014; Carvajal et al., 2016). 

Our morphometric research showed significant 
differences in forewing, 


hindwing, head and 
scutellum between the "rounded humerus form" and 
“ needle-like 
However, as that has been found for Anotogaster 
sieboldii (Kiyoshi and Hikida, 2012) , their centroid 


sizes do not differ greatly. Hence, an analysis of 


humerus form” of E. aeneus. 


shape may be more sensitive for detecting differences 
than that of size, because even when the size 
overlaps, the lineages differentiate. 

The morphometric patterns revealed in our study 
were discordant to previous molecular studies ( Zhao 
2015 ). While there were 
differences in morphology between the 
humerus form” and “needle-like humerus form” of 


et al., significant 


“ rounded 


the genetic differentiation on the 
mitochondrial sequence of the two forms was not 
distinct. Similar results were also found in Large 
Blue ( Maculinea arion) ( Bereczki et al., 2013). 
Mahalanobis and  Procrustes 
between the two forms for all the four characters were 
smaller than those interspecific ones ( Table 3). The 
DFA analysis gave similar results ( Table 4). Shape 
divergence showed by our study may not reflect 
specific divergence. 


E. aeneus, 


distance distance 


The origins of these differences may be resulted 
from some adaptive constrains such as the altitude, 
host plant, temperature, humidity, etc. ( Gumiel et 


al., 2003). It suggests that the size changes could 
be more adaptive than the shape itself ( Aytekin et 
al., 2007; Talloen et al., 2009; Outomuro and 
Johasson, 2011; Sanzana et al., 2013; Bai et al., 
2015). Because of the absence of allometry , and the 
overlapped centroid size, we supposed that the shape 
differences we experienced here should not be 
determined by adaptive factors. 

De Queiroz (2007) stated that species are 
equated with separately evolving metapopulation 
lineages, and many of the properties, such as 
morphological, genetic or behavioral differences, are 
more appropriately interpreted as operational criteria 
(lines of evidence) relevant to assessing lineage 
separation. The closer the lineages get to complete 
differentiation, the easier the diagnosis becomes 
( Wei et al., 2010; Marchiori et al., 2014). Hence, 
E. aeneus may be in the early stages of speciation. 

The two forms of E. aeneus can be separated 
clearly by all of the four characters in CVA with P- 
values for Procrustes distance and Mahalanobis 
distance all less than 0. 01. However, the correctly 
assigned percentages of the two forms in DFA were 
varying ( Table 4). This may indicate that the values 
of the four detect 
difference were not equivalent. The scutellum is the 
most informative equal 
intra- and 


characters to intraspecific 


character with almost 
correctly assigned percentages in 
interspecific level. 

It is a fact that molecular data always tend to 
receive more emphasis than morphological ones. 
However, not all of them are suitable for 
discrimination of all populations or subspecies 
( Sheppard et al., 1996; Bereczki et al., 2013; 
Oleksa and Tofiski, 2014; Robins et al., 2014; 
2018 ). 


morphometrics have clearly indicated its importance 


Hulme-Beaman et al., Geometric 
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in clarifying taxonomic and evolutionary relationships 
among taxa (De la Riva et al., 2001; Gumiel et al., 
2003 ; Campos et al., 2011; Bai et al., 2011, 2012, 
2014; Dehon et al., 2014) , and it was proved to be 
more useful than molecular ones in some studies 
(Bloch et al., 2010; Bereczki et al., 2013; Oleksa 
and Tofiski, 2014). In the present work, we are 
able to show that geometric morphometrics is 
efficient in describing intraspecific shape difference 
of E. aeneus. 
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摘要 :【 目 的 】 北 二 星 晴 Eysarcoris aeneus 是 一 种 广泛 分 布 于 古 北 区 的 重要 农业 害虫 ,可 危害 多 种 经 
济 作物 。 其 刺 肩 型 和 钝 肩 型 前 胸 背 板 后 侧 角 长 度 存 在 差异 ,暗示 其 种 内 变异 的 存在 。【 方 法 】 基 于 
几何 形态 学 的 多 元 回归 分 析 (multivariate regression )、 主 成 分 分 析 (principal component analysis, 
PCA), 32 # & & Ap (canonical variate analysis, CVA) 和 判别 函数 分 析 (discriminant function 
analysis, DFA) ,对 采集 于 中 国 19 个 地 区 的 142 头 标本 (98 3k 36 — X 3 AA Fe 44 头 作为 外 群 的 广 
— X Wi. E. ventralis 标本 ) 的 前 翅 、 后 翅 、 头 和 小 盾 片 4 个 性 状 进行 分 析 比 较 ,研究 北 二 星 晴 钝 户型 和 
刺 肩 型 之 间 的 形态 变异 。 【结果 】 对 于 所 研究 的 4 个 性 状 来 说 , 北 二 星 晴 钝 户型 和 刺 肩 型 标本 间 均 
未 检测 到 异 速 生长 的 存在 。 主 成 分 分 析 中 , 钝 肩 型 和 和 刺 肩 型 的 标本 均 有 重 壹 ,典型 变量 分 析 则 显示 
它们 之 间 存 在 显著 差异 ( 马 氏 距 离 和 普 氏 距离 的 已 值 均 小 于 0.01)。 判 别 函 数 分 析 显 示 , 依 据 这 4 
个 性 状 两 型 样本 间 的 正确 判别 率 在 67% ~89% 之 间 。【 结 论 】 结果 表明 ,前 翅 \ 后 翅 、 头 和 小 盾 片 都 
可 以 作为 北 二 星 暑 钝 户型 和 刺 户 型 之 间 形 态 变 蜡 的 评价 指标 ,小 盾 片 特征 具有 最 高 水 平 的 鉴别 价 
值 ;对 于 所 研究 的 这 4 个 性 状 来 说 , 北 二 星 晴 钝 户型 和 刺 肩 型 之 间 都 表现 出 显著 的 形状 差异 ,而 它 
们 的 质心 距离 差异 不 显著 ,说明 相对 于 体型 大 小 差异 分 析 ,形状 差异 分 析 能 更 灵敏 地 揭示 谱系 间 的 
变异 。 本 研究 表明 几何 形态 学 可 以 有 效 地 描述 北 二 星 晴 的 种 内 形态 变异 ,为 晴 类 昆 吕 种 内 变异 的 
研究 英 定 了 基础 。 

关键 词 : 北 二 星 晴 ; 几何 形态 学 ; 主 成 分 分 析 ; 典型 变量 分 析 ; 判别 函数 分 析 
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